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Recent direct numerical simulations of the FENE-P model of non-Newtonian hydrodynamics 
revealed that the phenomenon of drag reduction by polymer additives exists (albeit in reduced 
form) also in homogeneous turbulence. We introduce here a simple shell model for homogeneous 
viscoelastic flows that recaptures the essential observations of the full simulations. The simplicity of 
the shell model allows us to offer a transparent explanation of the main observations. It is shown that 
the mechanism for drag reduction operates mainly on the large scales. Understanding the mechanism 
allows us to predict how the amount of drag reduction depends of the various parameters in the 
model. The main conclusion is that drag reduction is not a universal phenomenon, it peaks in a 
window of parameters like Reynolds number and the relaxation rate of the polymer. 



PACS numbers: 47.27-i, 47.27.Nz, 47.27.Ak 

I. INTRODUCTION 

The phenomenon of drag reduction by polymer addi- 
tives is usually studied in channels or pipes, where the 
boundary conditions and the effects of the walls are very 
important til |4] • Until recently it was not known 
whether drag reduction can be achieved also in homoge- 
neous flows; this question has been answered recently in 
the afhrmative, via Direct Numerical Simulations (DNS) 
of the FENE-P model equations 0, |^ in homogeneous 
conditions (i.e. in a box with periodic boundary condi- 
tions) The FENE-P model takes the effect of the 
polymers on the Newtonian fluid into account by intro- 
ducing the conformation tensor R of the polymers into 
the fluid stress tensor. The FENE-P equations are known 
to model well the effects of polymers on the hydrody- 
namic flows, and DNS of these equations in channel ge- 
ometry recaptured very well the characteristics of drag 
reduction in experimental channel turbulence . The 
observation of drag reduction in homogeneous conditions 
offers an opportunity to investigate the phenomenon in- 
dependently of boundary layers and wall effects. Never- 
theless the FENE-P equations are relatively cumbersome 
to analyze without the help of DNS. The aim of this paper 
is to introduce a shell model of the homogeneous FENE- 
P equations. We will demonstrate that the shell model 
recaptures the main findings of the homogeneous DNS, 
and that these findings are understandable analytically, 
taking advantage of the relative simplicity of the shell 
model. To derive the shell model for drag reduction we 
make use of a formal analogy between the FENE-P equa- 
tions for viscoelastic flows and magneto-hydrodynamics 
(MHD) [13. It had been pointed out that if we form 
a tensor Ri^j from the direct product of the magnetic 
field Bi, i.e. Ri j = BiBj, then the nonlinear couplings 
of MHD lead to equations for the tensor R whose non- 



linear terms are equivalent to those of FENE-P, up to 
terms that remove the dynamo effect. This analogy is 
revisited and exploited in Sect. ^ The shell model for 
viscoelastic fiow is introduced and discussed in Sect. IIIII 
In Sect. II VI we present numerical simulations of the shell 
model and demonstrate the existence of drag reduction. 
In Sect. 0we present the mechanism of drag reduction. 
This is the central section of this paper. We show that 
drag reduction is not a universal phenomenon. Rather, 
it depends on the parameters, like the Reynolds number 
and the relaxation time of the polymer. The amount of 
drag reduction peaks in a window of these parameters. In 
Sect. VI we demonstrate that understanding the mecha- 
nism provides us with predictive power that we can test 
against numerical simulations. We conclude in Sect. IVIII 
by observing that precisely because drag reduction is not 
a universal phenomenon it can be manipulated by opti- 
mizing parameters. 



II. THE FENE-P EQUATIONS AND THEIR 
RELATION TO MHD 

The addition of a dilute polymer to a Newtonian fluid 
gives rise to an extra stress tensor T{r,t) which affects 
the Navier-Stokes equations j^] : 

— + {u-V)u = -Vp + i/sV^n + V-T , 

V • M = . (1) 

Here u{r,t) is the solenoidal velocity held, p{r,t) is the 
pressure and Vs is the viscosity of the neat fluid. In the 
FENE-P model the additional stress tensor T is deter- 
mined by the "polymer conformation tensor" R accord- 
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ing to 



T 



fir,t) 



pI 



R{r,t) - 1 



(2) 



Here 1 is the unit tensor, Vp is a viscosity parameter, r 
is a relaxation time for the polymer conformation ten- 
sor and Po is a parameter which in the derivation of the 
model stands for the rms extension of the polymers in 
equilibrium. The function /(r, i) limits the growth of 
the trace of i2 to a maximum value pm '■ 
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pI 



Pi 
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The model is closed by the equation of motion for the 
conformation tensor which reads 



dt 



[u ■ V)R^f 



R. 
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[f{r,t)R 
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■7/3 



(4) 



This model was simulated by DNS in channel flow tur- 
bulence, showing qualitative and quantitative agreement 
with laboratory experiments on drag reduction. Recently 
the same model has been used to understand whether 
or not drag reduction is observed in homogeneous and 
isotropic conditions 0- In homogeneous and isotropic 
turbulence, drag reduction can be determined by com- 
puting the ratio 



D 



cL 



(5) 



where E is the kinetic energy,e is the total rate of energy 
dissipation and L is the scale of the external forcing. The 
above expression of drag reduction can be easily reduced 
to the so called skin friction factor for turbulent chan- 
nel flows. The numerical simulation of homogeneous and 
isotropic turbulence were performed in a cube with pe- 
riodic boundary conditions. The external forcing was 
applied with random phase in order to ensure isotropy 
and homogeneity. The numerical simulations were per- 
formed for the Navier-Stokes equations and the FENE-P 
model for the same external forcing. Both the total en- 
ergy dissipation and the kinetic energy increased for the 
FENE-P as compared to the Newtonian case. A direct 
computation of D shows that there is a drag reduction of 
about 20%, i.e. roughly of the same order as what had 
been observed in turbulent channel flow. Also, in homo- 
geneous and isotropic turbulence, the Taylor microscale 
appeared to increase, apparently precisely as much as 
the buffer layers increases in channel flows 0]. This is 
an interesting result because it tells us that the effect of 
boundary conditions is not crucial for drag reduction, at 
least from a physical point of view. Nevertheless, it is 
still difficult to understand from numerical simulations, 
even in the homogeneous and isotropic case, what is the 
physical mechanism that is responsible for drag reduc- 
tion. The increase of the Taylor microscale is certainly 



not enough to explain quantitatively the increase of the 
kinetic en ergy, as somehow previously suggested in the 
literature (JjlJ. 

Having understood that the homogeneous simulations 
exhibit drag reduction, we would like to propose a mecha- 
nism for it. Rather than doing it directly with the FENE- 
P model, we would present first a simplified model. We 
have already shown before that drag reduction appears 
in simplified models like the Burgers equation . Here 
we derive a shell model for the FENE-P equations. The 
advantage of the shell model is that it is much more 
tractable analytically than the full FENE-P equations. 
We will present the model, demonstrate explicitly that 
it exhibits drag reduction in much the same way as the 
FENE-P equations, and finally offer a new mechanism to 
understand the phenomenon. 



III. THE SHELL MODEL 

To derive a shell model of the homogeneous FENE-P 
equations (without boundaries) we proceed in two steps. 
First we recall a recent remark jl^ that the FENE-P 
equations can be recaptured almost entirely by taking 
the conformation tensor to be a diadic direct product of 



of a vector B, i.e Rij 
the equations read 



BiBj. In terms of this vector 



du 
'dt 



+ {u ■ V)u 



-Vp+B ■ (VB) + UsV^u 



V u = , 



— + B (Vm) 

T 



V B = 



(6) 



These equations are identical to the FENE-P model up 
to the explicit appearance of the function f{r,t). The 
learned reader of course recognizes that for r — > oo these 
equations are isomorphous to magneto-hydrodynamics 
(MHD). We can therefore write immediately, by inspec- 
tion, a shell model for FENE-P by using the well studied 
shell model for MHD [isL , including the relaxation 
term for finite r. We denote the velocity field by w, 
and the "polymer " field by B. The dynamical variable 
of the shell model are the field at wave vector fc„, de- 
noted respectively as w„ = and i?„ = i?(fc„). The 
shell model restricts attention to wavevectors A:„ = fcgA", 
where typically in numerical simulations A = 2. 

In order to derive the shell model equation, we consider 
the following non linear operator: 

$„(m, B) = 6ifc„W„+2S* + l - &lfc„-lMn+l-B*_i -I- 
Clfc„M*_|_ii?„+2 + Cikn-2Un-lBn-2 + 

2Un-2Bn-l (7) 

where fei = 1 — 6; ci = 1 + c; 5c = 6 + c, and —1 < 6 < 
and c — \ + h are the usual parameters defined in the 
Sabra model. 
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In terms of the non linear operator $, the Sabra shell 
model of turbulence |115| can be written as 

= g^'ri(w, U) - vk^Ur, + /„ (8) 

where /„ is an external forcing and v the kinematic vis- 
cosity of the model. Let us remark that the following 
relation can be proved: 



il]„$„(u, B)Bl - il]„$:(u, B)Br, 



0. 



(9) 



Using the non linear operator $ it is possible to model 
equations ® in the framework of shell models, namely: 



^$„(m, u) - ^$„(B, B) - Uklu + fn 



dB„ 



Mu,B) - -$„(B,u) 



1 



-B„ 



(10) 



Equation @ tells us that the generalized energy E, 

E = Eu + Eb 

Eu = ^nUnK - Eb = ^nBnB* . (11) 

oo and 



is conserved in the inviscid limit, i.e. for r 
1/ ^ 0. 

We will refer to this model as the SabraP model. Be- 
side the generalized energy E, the model conserves the 
"cross helicity" in the inviscid limit 



{u„Bn) 



(12) 



In MHD one needs to worry about the existence of a 
dynamo effect, i.e. an unbounded increase in the mag- 
netic field. In our case the term that models the polymer 
relaxation time —B^/t will be responsible for guarantee- 
ing stationary statistics without dynamo. In addition to 
the conservation laws the equations of motion remain in- 
variant to the phase transformations w„ —>■ Un exp(i0„) 
and Bn — + Bnexp^itpn). The conditions are 
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Tpn - 


\- 4'n+l 


- <i>n+2 


= . 


(16) 



This implies ipn — 4'n Vrt. As a result of the phase con- 
straints there exist in this model only few non-zero cor- 
relation functions. The only second order quantities are 
and The only third order quanitites are 

of the form 3(/3„_i/3„/3*_^j) where /? can be u or B. 



IV. NUMERICAL INVESTIGATION OF THE 
SHELL MODEL: DRAG REDUCTION 

In this section we compare the solutions of the shell 
model ((Tn|) to the usual Sabra shell model for the cor- 
responding Newtonian flow. The Sabra model ijHJ) is ob- 
tained from H1U|) in the limit r ^ 0. Alternatively, we 
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FIG. 1: Kinetic energy of the SabraP model for v = 10~^ 
as a function of r. The constant reference line corresponds 
to the kinetic energy computed for the Sabra model without 
polymer. 



can get the Sabra dynamics by simply taking as initial 
conditions i?„ = 0. 

To have a meaningful comparison we always drive the 
two models with a constant power input. In other words, 
we choose 



« Fa p E\y 

J \ ^ J 2 J. 



fn ^0 for n>2 



(17) 



with Fa — Fb = 10"'^ (2 -I- 2i). Since the power input is 
the same, drag reduction is exhibited in if the kinetic 
energy of the flow increases. The latter is simply 
We will investigate the existence of drag reduction, its 
dependence on parameters, the question of the dissipative 
scale, and the dynamical signatures of drag reduction. 



A. Drag reduction and its dependence on 
parameters 

We have numerically investigated the behaviour of the 
SabraP model for different values of r and v. In figures 
(P)-© we show {Eu) for three values of the viscosity and 
for different values of r. For concreteness we have fixed 
the model parameter 6 to be —0.4 in all the simulations. 
In all the figures the constant line corresponds to the 
value of the kinetic energy computed for the Sabra model 
without coupling to Bn- By inspecting the three figures, 
one can safely state that the SabraP model shows drag 
reduction. In particular, for all cases, there is an optimal 
choice of r for which the effect of drag reduction is max- 
imal. For r — > and r — > oo drag reduction decreases 
and eventually we enter a region of parameters where we 
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FIG. 2: Kinetic energy of the SabraP model for v — 10~ 
as a function of r. The constant reference line corresponds 
to the kinetic energy computed for the Sabra model without 
polymer. 
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FIG. 3: Kinetic energy of the SabraP model with polymer for 
V — 10^* as a function of t. The constant reference line cor- 
responds to the kinetic energy computed for the Sabra model 
without polymer. 



B. Which scales are responsible for drag reduction? 

To understand which scales are responsible for the 
drag reduction, we compare (|w„p) for both models with 
V = 10~^, again at the same power input. This is shown 
in Fig. ^ This figure teaches us an interesting and im- 
portant lesson. It is clear that the drag reduction is due 
to the relative increase in (IwriP) /or small values of n, 
and that this occurs on the expense of a relative decrease 
in (|u„p) for high values of n. This finding is in close 
correspondence with similar conclusions obtained for the 
FENE-P model, both in homogeneous and channel flows 



C. The dissipative scale 

In some theories of drag reduction it was proposed that 
the dissipative scale is increased in the viscoelastic flow, 
and that somehow this is responsible for the phenomenon 
I1II1II3- To test this possibility we plot in Fig. IHthe 
quantity as function of n. This quantity peaks 

at the dissipative scale, i.e. the Kolmogorov scale. In- 
specting figure (jSJ teaches us that the dissipative scale 
has not changed at all between the Sabra and the SabraP 
models, even though the latter certainly exhibits drag re- 
duction. Thus, as indicated before, drag reduction should 
be understood as a phenomenon of the energy containing 
scales rather than the dissipative scales. 



observe drag enhancement. Moreover, for fixed value of 
T and decreasing values of drag reduction decreases, 
reaching a mere few per cents for v = 10~®. 



D. Dynamical signature of drag reduction 

The similarity between the FENE-P and its shell ana- 
log transcends statistical quantities. To observe the 
close dynamical similarity it is instructive to consider the 
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MECHANISM FOR DRAG REDUCTION 




FIG. 5: Energy dissipation computed for each shell for the 
Sabra model (dotted line with circles) and for the SabraP 
model (continuos line with squares) 




FIG. 6: Time behavior of the quantity 11 as defined in the 
text 



quantity 



n = zS„<$„(B, B) ~ zI]„?/„$;(B, B) , 



(18) 



which describes the exchange between the kinetic energy 
Eu and the "polymer" or "elastic" energy Eb- In figure 
© we show a time series of 11 for ly = 10^^. 11 is always 
negative; the effect of the "polymers" is to drain energy 
form the kinetic energy. Moreover the dynamics of 11 is 
strongly intermittent which is a feature already observed 
in the DNS of the FENE-P model. The numerical simu- 
lations indicate the conclusion that the model introduced 
in this paper shows drag reduction in a way qualitatively 
close to the observed behaviour of the FENE-P model 
i 

We note in passing that it is not the first time that shell 
models seem to reproduce many of the features of turbu- 
lent flows; it is gratifying however that we can present a 
similar success even when we include relatively non triv- 
ial effects induced by polymer dynamics. 



This is the central section of this paper, in which we 
propose a detailed mechanism for drag reduction in the 
present model. We begin by analyzing the necessary con- 
ditions for drag reduction. 



A. Necessary condition for drag reduction 

To derive a necessary condition for drag reduction, let 
us consider the equation for the total energy, which reads: 



1 



dE/dt = S„ ^ (/„< + f*u„) - vklunU*^ 

BnB* 



(19) 



At steady state, with power input maintained constant 
at P, we have 



T 



(20) 



All the terms on the RHS are strictly positive. Since 
the energy input P is constant for the Sabra and the SP 
models, we get 



^nkl [(u„Os - (w„Osp] > 0. 



(21) 



On the other hand if the SP model is to be drag reducing, 
we must have 



S„ [(w„u* )s - (w„u* )sp] < 0. 



(22) 



The only way 121|) and H22|) can hold simultaneously is 
if for small fc„, |w„|5p > \un\s, and sufficiently larger to 
compensate for the fact that at large fc„, |w„|sp < junls. 
This means that the kinetic energy plotted versus k has 
to display an increased slope at least somewhere for drag 
reduction to take place. We have seen this already in Fig. 
21 We show this important phenomenon once more in a 
log-log plot in Fig. [7| in which also the i?„-spectrum is 
shown for future reference. We see very clearly the cross- 
ing that occurs between the Un spectrum of the SabraP 
model and the Sabra counterpart, which is the necessary 
condition for drag reduction. Note that the increase in 
slope is a necessary but not a sufRcient condition for drag 
reduction. We may increase the slope but not enough to 
cross the Sabra spectrum, or cross but not far enough to 
compensate for the reduced kinetic energy at large k. 



B. Typical scales related to the polymer 

A discussion of the mechanism of drag reduction calls 
for pointing out the existence of two typical scales that 
were already introduced in the past in the literature on 
drag reduction. The first is the Lumley scale, kc, which 
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FIG. 7: Energy spectrum of the SabraP model (continuos 
line with squares) and the Sabra model (continuos line with 
circles) for v — 10~^. The continuos line with black triangles 
represents the energy spectrum of the B field. 



FIG. 8: The inverse of the "eddy turn over time" ^ {En)kn as 
a function of k„ (continuous line with squares) . The constant 
reference dotted line is r^^. The crossing in inertial range 
identified kc- In this figure u — 10~^, r — 0.4. 



is defined by the relaxation time of the polymer being 
of the same order as the eddy turn over time. For our 
model this scale satisfies 

u{k,)k, - r-i . (23) 

Note that by definition this scale is Reynolds number 
independent. 

The other scale, that we refer to as the de Gennes scale 
kg, is where the kinetic energy on the scale kg is of the 
same order as the elastic energy: 

U^{kg)r^B^kg) . (24) 

In fact, in the SabraP model the scales so defined appear 
to be very close, if not identical to each other. In par- 
ticular, we will show presently that also kg is Reynolds 
number independent. To demonstrate the equivalence of 
the two scales we first exhibit in Fig. IHlthe numerical es- 
timate of kc- The physical significance of kc is not in the 
accidental identity of two time scales, but rather that for 
A;- vectors smaller than kc the effect of the Bn field on the 
energy flux is negligible, but not so for /c- vectors larger 
than kc- To see this introduce two quantities related to 
the energy flux in the SabraP model, namely: 

Sn = (/to«_iU>„+i)) (25) 
r„ = (/m(i3:_i<B„+i)) (26) 

The physical meaning of the two quantities is rather 
clear: S'„ describes the flux of kinetic energy from large 
scale to small scales due to non linear terms, while r„ 
describes the flux of kinetic energy to the polymer field. 



We expect that for fc„ near k^ the effect of T„ cannot be 
neglected in the dynamics, i.e. the average energy fiux 
for the velocity field Gn = Sn — Tn begins to change with 
respect to what it is observed in the Sabra model. 

In figure © we show the quantity Gn computed for the 
same model parameters of figure ^ . The symbols refer 
to the Sabra model while the continuous line corresponds 
to the SabraP model. In the vicinity of Uc 5.5 the two 
models show a different behaviour and in particular the 
SabraP model shows a decrease of the total energy flux 
Gn as previously claimed. 

Regarding the scale kg, it can be read from the spec- 
trum shown in Fig. |7| in which = 10~^. In Fig. 1101 we 
show the analogous spectra for ly = 10^^. Clearly kg did 
not change at all, in agreement with our assertion that it 
is Reynolds independent. Finally, we note that in all the 
flgures shown kc and kg are of the same order of mag- 
nitude, and in the sequel we do not distinguish between 
the two. 



C. The effect of the polymer at large fc-vectors 

In this subsection and the next we discuss the effect of 
the En fleld on the u„ field for fc- vectors much larger and 
much smaller than kc- We will show that the spectrum 
(|u„P) exhibits essentially the same scaling exponent as 
the Sabra model, but the amplitude is affected by the 
presence of the Bn field. This will be an important in- 
gredient in the mechanism of drag reduction. 

Begin with fc„ large, fc„ 3> fcc. In this regime the effect 
of the relaxation time r on the dynamics of the Bn field is 
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FIG. 9: The average energy flux Gn computed for the SabraP 
model (continuous line with squares) and the Sabra model 
(dotted lines with circles). 
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FIG. 10: Energy spectrum of the SabraP model (continuous 
line with squares). The line with black triangles represents 
the energy spectrum of the B field. 



ilar slopes, although intermittency affects the two spectra 
in a different way. 

On the other hand, the amplitudes of the two spec- 
tra need not be the same. The relative displacement of 
the two power laws is determined by numerical details in 
the model. To estimate this displacement we will esti- 
mate the amplitudes of the two spectra at the dissipative 
scale. The contribution to the dissipation of u is mainly 
from the small scales, i.e. very large values of kn- We 
can define an effective scale kd the scale at which energy 
dissipation peaks: 

mldunf) ^ kj{\udf) (27) 

where Ud = u{kd), and kd is of the order of the Kol- 
mogorov scale. Since we found that the dissipative scale 
hardly changes when we add the coupling to the i3„ field, 
we can deduce from (|2U|I that 

{\udf)sP-{\udf)s^T.^^- (28) 

n 

On the other hand, the sum on the RHS of Eq. 128|l is 
a geometric sum dominated by the contribution of B{kc) 
where Bn is maximal. We thus estimate the relative dis- 
placement of the two spectra at high values of fc„ by 

{\ud\')sP^{\ud\')s^^^^^. (29) 

Thus to first approximation we expect the slopes of the 
two spectra to remain unchanged, maintained at a con- 
stant difference from each other as given by (|29|l . until fc„ 
approaches k^ from above, where the effect of the relax- 
ation time T on the dynamics of cannot be neglected. 

D. The effect of the polymer at small fc- vectors 

Next we discuss the slope of the u„ spectrum for /c„ ^ 
kc- This is very easy, since the amplitude of Bn is very 
small due to the very efficient exponential damping by 
T. Thus the w„ field hardly feels the coupling to Bn, 
and its slope, up to intermittency corrections, is again of 
the order of C2 = 2/3. Again, the amplitude is changed 
compared to the pure Sabra case, and this is the most 
important feature that is discussed next. 

E. The tilt in the spectrum at /c„ « kc 



completely negligible. The dynamics of _B„ is dominated 
by its coupling to u„, simply because Unkn > T"^. But 
then in this regime the dynamics is like the one of MHD 
which had been analyzed in detail in The central 

conclusion of that analysis is that up to intermittency 
corrections, the spectra of both the Bn and the m„ fields 
exhibit a scaling exponent C2 = 2/3. Indeed, inspecting 
Fig. ^1 we see that for large /c„ the two spectra have sim- 



Considering the spectrum in Fig. IIUI we note that the 
Bn spectrum increases rapidly when kn kc from the 
left. To understand this phenomenon consider the equa- 
tion of motion for i?„ at steady state. To leading order 

= (1^1) - (|cifc„<+iS„+2|) - ^ , (30) 

where we have neglected terms of the order of i?„+i, but 
including them will lead to similar conclusions. Using the 
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fact that |fc„+iM*^_;^T| <C 1, and since |fc„+iu*+;^i3„+2| < 
|fc„+iM*_|_]^||-B„+2| we immediately conclude that 



|B„|) « (|B„+2|) 



(31) 



We continue this argument recursively to estimate the 
largest polymer contribution B{kc) as 



{\Bo\) 



{\hui\){\k3U3\)...{\ke-lU,-i\)T^'/^ 



(32) 



where A""^ = fee- 

In the vicinity of the scale kc, we have to leading order 
in B in the kinetic energy equation, 



= -fc„5„+i - bknSn + (1 + b)kn^2Sn- 



(33) 



When the amplitude of the polymer goes to zero {Bn — > 
Vn) the only solution is the well known scaling law 
Sn oc k^^. However the last term in H33I) forces now a 
tilt in the spectrum. Its sign is exactly such that Sn-i 
has to increase compared to Sn and respectively S'„+i. 
Of course, for fc„ <C kc the effect of the i?„ field on the 
M-spectrum is again negligible, and therefore the spectral 
slope will settle back to the Sabra value. However if 
the tilt in the vicinity of kc results in crossing the Sabra 
spectrum we would have a whole spectral range where 
the energy is higher. 

We therefore conclude that the existence of drag re- 
duction depends rather heavily on the sign of the energy 
transfer at scales close to kc. To check the sign directly 
in the numerics and thus to substantiate the existence of 
the tilt we return to the equations of motion and write 



, ^„ , (u, u, u) - {u, B,B)- vki\u„ 

|2_,T,(")/„, u U\ I ,T/(")^„. D ON -'-ID |2 



-,(«) 



dt 



\B„V = {u, B,B) + (", B,B)--\B„\' , (34) 



where the term ^I* J"-* (u, w, u) represents the kinetic energy 
flux of the field, while '^'^\u, B, B) is the energy flux 
going from the velocity field to the polymer field. Finally, 
term "i'^^\u, B, B) is the flux of energy of the polymer 
field due to the transport of the velocity field. Figure 1)11(1 



shows ^ 



: 1 , ^ ^"^ and 4-^"^ for r = 0.4 and 1/ = 10^ the 
same parameters of figure {T)). It is important to observe 

(n) 

that '^l becomes positive for n > Uc- For a given n the 

term can be written as 5'^"'' — Ln~ Sn where L„ is 
the amount of energy flux given from the large scale to 
scale kn and Sn is the amount of energy flux given from 
scale kn to smaller scales. It follows that when the energy 
flux is constant L„ — Sn and therefore \E'^"''(n) = 0. On 

the other hand, a positive value of ^'^"■'(ri) implies that 
Ln > Sn- This is exactly what is shown in flgure (|ll|l . 
The imbalance of the energy flux is compensated 
by the flux of energy from u„ to Bm given by the term 

in) 

—^2 ■ It interesting to observe that the last term in 




FIG. 11: Energy fluxes for the SabraP model. The continuous 



line corresponds to the fluxes of the kinetic energy The 



squares correspond to — ^ 
triangles is ^fg"' 



The dashed line with black 



(n) 

the balance equation, namely , is rather small, i.e. 
the effect of an energy cascade of the polymer is rather 
weak. 



F. Discussion 

While we have been able so far to describe a convinc- 
ing scenario for drag reduction, we still should explain 
the mechanism for the increase of the large scale energy. 
Since the field _B„ is negligible for small n, the average 
energy fiux per unit time at small n must equal the in- 
put work per unit time at the largest scales. However, the 
energy flux does show time and scale fluctuations which 
could behave differently for the Sabra and SabraP models. 
More specifically let us consider the quantity G„ defined 
in Subsect. IV Bl As already discussed, G„ represents the 
energy flux at scale fc„ due to both the non linear terms 
in the velocity field and the non linear term in the Bn 
field. In terms of G„, we can build a large scale energy 
flux Wl = G2 + G3 + G4 which represents the full amount 
of energy flux across the largest scales, namely across k- 
vectors kn < kc, for which the average energy flux is 
invariant to changing Sabra to SabraP. The definition of 
Wl is such that Wl > means an energy flux from large 
scales to small scales. In figure (I12|l . we show the proba- 
bility distribution of for both models, with numerical 
parameters = 10~^ and t = 0.4. The vertical line in the 
figure indicates the average value, which, as expected is 
invariant. As one can observe, the two probability distri- 
butions show a substantial difference for negative values 
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-0.02 0.02 0.04 0.06 900 950 1000 1050 1100 



n L 

FIG. 12: Probability distribution functions of Wl for Sabra FIG. 13: Time series of the kinetic energy (continuous line) 

model (line with black triangles) and SabraP model (contin- and 2 x Wl (dotted line) for the Sabra model, 
uous line). 



of Wl- in the SabraP model one can have larger and 
more frequent negative values of the energy flux. This 
can happen even if the instantaneous value of the total 
energy per unit time obtained by the polymer field from 
the velocity field is always positive! Qualitatively this 
means that while the velocity field is always forcing the 
polymer field, at very large scales the flux can, from time 
to time, be reversed, and the polymer field forces there 
the velocity. This is how the amplitude of the energy 
spectrum is being increased on the average. 

Clearly, this mechanism could not work unless the 
"forcing" by the polymer field acted in phase with the 
growing kinetic energy. In order to clarify this further 
we present two time-series of the kinetic energy and Wl , 
in Fig. H13|) for the Sabra model, and in Fig. (|14|) for the 
SabraP model. A close inspection of the figures shows 
that the reverse of the energy fiux Wl occurs exactly 
during the growing phase of the kinetic energy, leading 
therefore to a larger value of the instantaneous kinetic 
energy. This in phase mechanism is responsible for drag 
reduction. Note that this mechanism strongly depends 
on the large scale dynamics and the value of kc- For fc„ 
larger than kc no significant difference in the statistical 
behaviour of the energy flux is observed. However, the 
amount of energy forcing, due to the polymer at large 
scale, can depend on the Reynolds number, at least in 
the SabraP model. If this is the case, then drag reduc- 
tion should depend on the Reynolds number only through 
the two relevant scales appearing in the systems, namely 
kc and Xt, the latter being the Taylor microscale 



Eu 




-0.025 - 



1100 



(35) 



FIG. 14: Time series of the kinetic energy (continuous line) 
and 2 x Wl (dotted line) for the SabraP model. 



which also depends on r in the SabraP model. Because 
the drag is a dimensionless quantity, we argue that the 
only way in which the Reynolds number may appear in 
the drag reduction is by means of the dimensionless quan- 
tity ij, = kcXriT). 

There is a simple argument, proposed in the next sec- 
tion, which explains why at large Reynolds numbers one 
may observe the same qualitative mechanism, i.e. drag 
reduction, with smaller effects on the kinetic energy. As 
a matter of fact, the numerical results show that drag 
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FIG. 15: Schematic of the effect of polymers in drag reduc- 
tion on the turbulence energy spectrum. Dotted line: neat 
fluid. Solid line: polymeric solution. The spectral slope is 
unchanged for large and small scales, while at scale kc there 
is a significant upward tilt. 



reduction reaches its niaximum for ~ 1 



VI. PREDICTIONS OF THE THEORETICAL 
MECHANISM 

We summarize the mechanism of drag reduction using 
the cartoon shown in Fig. E| The tilt in the spectrum 
occurs in the vicinity of with the asymptotic slope for 
kn kc and fc„ 3> kc remaining essentially unchanged. 
With such a spectrum the two inequalities H21|) and H22|) 
are obviously obeyed. 

The difference in the spectra for fc„ ^ kc is deter- 
mined predominantly by Eq. H28(l . This equation pre- 
dicts that this difference will be greatly increased when 
the Reynolds number is increased (i.e. when z/ ^ 0), see 
Fig. El Of course, if this happens we can lose the whole 
effect of drag reduction, since the amount of tilt at kc 
is basically independent of v. We need to maintain the 
spectral difference small enough for the tilt to effect a 
crossing of the spectrum of SabraP and Sabra. Also the 
position of kc is important. If we reduce kc (i.e. increase 
r) the tilt is too far to the left and therefore it will fail 
to increase the energy. In fact it can be drag enhancing. 
The combined effect of decreasing v and increasing t is 
shown in Fig. El Needless to say, also if we decrease r 
too much we may lose drag reduction since the tilt will be 
pushed to the irrelevant dissipative range where no en- 
ergy containing modes exist. Also, if r becomes too low, 
Bn becomes smaller, and the amount of tilt is decreased, 
as can be seen directly from Eq. Although decreas- 

ing the field i?„ brings the spectra closer together in the 



FIG. 16: Schematic of the turbulence energy spectrum when 
the polymer relaxation time scale is too slow, and the 
Reynolds number is too large. 




K 



FIG. 17: The relaxation time r is to low. 



large fc„ regime, the tilt may not suffice to reduce the 
drag. Such a situation is shown schematically in Fig. El 
Actually, using the language introduced in the previous 
section and the above discussions, we are able to give an- 
other argument to understand how drag reduction could 
depend on the Reynolds number. As previously said, the 
relevant dimensionless number in the system is ^ = kcXr- 
If r — > C!0 then fee — > and we know that drag reduction 
must be inhibited. It follows that for /i ^ we can- 
not observe drag reduction. For fixed kc and increasing 
Reynolds, Xt{t) decreases as well, although not necessar- 
ily as i?e~^/^, where Re is the Reynolds number. Then, 
for fixed r and increasing Reynolds number we should ob- 



11 



serve a decreasing effect of drag reduction, as observed 
in our numerical simulation. The same reasoning can be 
applied to get information for small Reynolds numbers, 
as the following argument shows. For t — > we have 
already shown that no drag reduction is possible simply 
because fcc oo. This is equivalent to say that when /i 
becomes too large there cannot be any drag reduction. It 
follows that for small Re, i.e. for large At, drag reduction 
disappears. 

VII. CONCLUSIONS AND DISCUSSIONS 

In this paper we discussed several points concerning 
the possible formulation of a theory for drag reduction 
in turbulent flow with dilute polymer. It is worthwhile, 
therefore, to review the main points. 

A) We have introduced a shell model resembling the 
dynamical properties of the FENE-P equations. Beside 
any theoretical considerations, the model shows drag re- 
duction in a way close to what already observed in the 
numerical simulations of the FENE-P model. The impli- 
cations of this result is that one need not focus on bound- 
ary effects or dynamical properties of coherent structure 
in order to capture the basic physics of drag reduction. 

B) There exist a relevant scale in the system, fcj, defined 
by the so called 'time criterion', i.e. u{kc)kc ~ . In 
the vicinity of this scale there is a tilt in the spectrum 
which causes a crossing of the SabraP velocity spectrum 
above the Sabra spectrum for fc„ < kc- This in its turn 
means an increase in the kinetic energy at large scales. 



Drag reduction can be physically understood in terms 
of the energy exchanges between the velocity field and 
the polymer field for fc„ ~ kc- We have succeeded in 
proposing a coherent picture, based on the equation of 
motions, for the dynamics which is in close agreement 
with the numerical results. 

C) Drag reduction is a property of large scale flow and 
its dynamics. This implies that a quantitative descrip- 
tion of drag reduction must depend on the details of the 
flow, the forcing mechanism as well as the Reynolds num- 
bers. Although a general qualitative mechanism should 
occur in all drag reduction flow, the amount of the drag 
reduction itself depends on how much energy is inter- 
mittently given to large scale velocity. Thus large scale 
fluctuations are important for a quantitative theory. 

D) Drag reduction by no means could be reduced to 
the dynamics at the dissipation scale. Although drag 
reduction could be Reynolds dependent, drag reduction 
cannot be reduced to a simple increase of the dissipation 
length. Actually, the dissipation scale does not seem to 
be affected by drag reduction. 
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